X-wave mediated instability of plane waves in Kerr media 
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Plane waves in Kerr media spontaneously generate paraxial X-waves (i.e. non-dispersive and 
non-diffractive pulsed beams) that get amplified along propagation. This effect can be considered 
a form of conical emission (i.e. spatio-temporal modulational instability), and can be used as a key 
for the interpretation of the out of axis energy emission in the splitting process of focused pulses in 
normally dispersive materials. A new class of spatio-temporal localized wave patterns is identified. 
X-waves instability, and nonlinear X-waves, are also expected in periodical Bose condensed gases. 



The dynamics of focused femtoseconds pulses (FTP) 
in optically nonlinear media has a fundamental impor- 
tance and it is relevant in all of the applications of ul- 
trafast optics. The basic mechanism, when dealing with 
the propagation in normally dispersive materials, is the 
splitting of the pulse, which has been originally predicted 
more than ten years ago P, |3| and recently reconsid- 
ered, due to the develo pme nt of the physics of FFP. 

If z is the direction of propagation, t is the time in the 
reference frame where the pulse is still, and r = \Jx^ 
the radial cylindrical coordinate, this process can be 
roughly divided into a series of steps. They describe 
the propagation of a gaussian (in space and time) pulse, 
which travels in a focusing medium, and undergoes rel- 
evant reshaping, due to the interplay of diffraction, dis- 
persion and the Kerr effect (when the peak power is suffi- 
ciently greater than the critical power for self-focusing 
HQ): 1) the out of axis energy is focused towards r = 
around i = 0; 2) the pulse at r = is compressed; 3) 
lobes appear in the on-axis temporal spectrum; 4) the 
pulse splits in the time domain. 

Before the breakup, relevant out of axis energy emis- 
sion and re-distribution occurs, as originally described by 
Rothenberg in Ref. 0. Looking to the spatio-temporal 
profile an X-shape (or hyperbola) is observed before 
the splitting (see, e.g., figures in 0, This pro- 

cess has been theoretically described by different ap- 
proaches 0, 0, Q, 0|, and the onset of an X-shape 
can be ultimately related to the hyperbolic character- 
istics over which small perturbations evolve. 0| This 
can be checked, for example, by the hydrodynamical ap- 
proach to the nonlinear Schrodinger equation. 0, 0| 

Recently, attention has been devoted to the existence 
of X-waves in optically nonlinear media. The latter 
are self-trapped (i.e. non-diffracting and non-dispersive) 
waves that are well known in the field of linear propaga- 
tion in acustics 19, 20, 21. 22) and in electromagnetism 
m m il S m Slalol- The simplest X-wave 
has the shape of a double-cone, or clepsydra, that ap- 
pears as an X when, e.g., a section is plotted in the plane 
(x,t), and as V in the plane (r,t), as shown in figure ^ 
Optical X-waves in nonlinear media have been theoret- 



ically predicted in |31| and experimental results, with a 
direct observation of the conical spatio-temporal shape, 
have confirmed their existence and generation in crystals 
for second-harmonic generation (SH-G) [23)IUI3- 
Ref. [3^ it has been theoretically shown how, during non- 
depleted-pump SHG, the phase matched spatio-temporal 
harmonics let the SH beam become an X-wave. 

Conical emission (CE), or spatio-temporal modula- 
tional instability_(MI) (Hi, S HI El), has been ad- 
dressed in Ref. |39| as a basic mechanism underlaying the 
spontaneous formation of an X-wave and, as a founda- 
tion for the understaing of the splitting, in Ref. 10. In 
this Letter I consider the formation of X-waves in Kerr 
media (or in quadratic media, in the regime where they 
mimic cubic nonlinearity, see, e.g., the chapter after Tor- 
ruellas, Kivhsar and Stegeman in and Ref. IZlj) . A 
new form of instability of plane waves can be introduced 
by directly involving self-localized spatio-temporal wave 
patterns. The process strongly resembles CE, i.e. the 
amplification of plane waves from noise, but in this case 
X-waves, instead of periodical patterns, emerge from the 
breakup of an unstable pump beam. This mechanism 
is responsible of the first stage of the pulse-splitting, 
i.e. the out-of-axis energy redistribution, such as MI 
breaks a continuos wave signal into a periodical pattern 
of solitons. 42] 

The wave equation describing the propagation in non- 
linear Kerr media, whose refractive n index obeys the law 
n — riQ + n2l, with / the optical intensity, in the frame- 
work of the paraxial and of the slowly varying envelope 
approximation, is written as 



idzA + ik'drA + 7;^ ^xr 
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where A is normalized such that \A\'^ ~ I, and I have 
taken k" > (i.e. the medium is normally dispersive). 
X, Y, Z, T are the real world variables, A the wavelength, 
n(A) the refractive index and fc(A) = 27rn(A)/A. Eq. 
can be casted in the adimensional form: 



idzU + Aj_u — dttu + x|upM — 0, 



(2) 



by defining z — Z/Zdf with Zdf = 2fcWo and Wq a ref- 
erence beam waist; Aj^ = dxx + dyy the tranverse Lapla- 
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FIG. 1: Tri-dimensional plot of the simplest radial- 
symmetry X-wave (A = 1). 



cian with {X,Y) = Wo{x,y); t = [T - Z/Vg)/Ta the 
retarded time in the frame travelhng at the group ve- 
locity Vg = 1/k' in units of Tq = (k" Zdj jifl'^ . The 
optical field envelope A is given by A = Aqu, with 
^0 = (no/fc|n2|^d/)^/^- X > (X < 0) identifies a fo- 
cusing (de- focusing) medium being n2 > (jii < 0). 
I start considering paraxial linear (i.e. X = 0) X-wave 
solutions of (0) (paraxial X-waves have been considered 
in detail in Ref. ^3) , defined by 



O = (Ai - a«)^ - 0. 



(3) 



Introducing the complex variable v = {A — it)^ + r'^, with 
A a real valued arbitrary coefficient, I have from (PJ the 
equation 



(4) 



from which ip = Ci/ ^/w -I- C2 . Ci and C2 are arbitrary 
complex coefficients. Note that both the real and the 
imaginary parts of this solution, are real-valued X-wave 
profiles. The former being the simplest X-wave, given by 
(the branch cut for the square root is along the negative 
real axis) 
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Remarkably V'x still holds when referring to the Helmotz 
equation, instead of the paraxial wave equation (see e.g. 
Ref. l43). Its plot is given in figure ^ 

X-wave instability can be introduced in the same way 
as MI, i.e. by perturbing the plane- wave solution of jSJ: 
a ~ ap = aoexp(iaQz), with ag a real-valued constant. 
Letting a = (oq + e(a;, y, t, z)) exp(iagz) I have, at first 
order in e: 

z9,e + Oe + X«o(e* + e)- (6) 
Writing e — e(r, i, z) = ipxi^, i)/^(-2) gives 

i5,M-f Xao(M + M*) = 0. (7) 



The perturbation is thus 

e = [a + i{(3 + 2xa^az)]il)x{r, t) 



(8) 



with a and (3 arbitrary real-valued constants. Equa- 
tion ijSJ) represents an X-wave which grows linearly along 
propagation (independently on the sign of x)j with am- 
plification given by the intensity of the plane wave ap. 
This situation strongly resembles MI, where plane waves 
are exponentially amplified at the expense of the pump 
beam, with a gain determined by the pump intensity. For 
this reason, it is natural to refer to this process as X-wave 
instability. As for MI, the amplified wave can be artifi- 
cially externally feeded, or it can be generated by noise. 
[igf Note also that X-wave instability can be triggered 
by conical emission. Indeed, as shown in [s^, the lat- 
ter generates the required spectrum to form an X-wave, 
which eventually gets amplified, as discussed above. 

The previous treatment can be generalized in several 
ways. In particular it is possible to show that exponential 
amplification of X-wave-like beams can be attained. It is 
necessary to enlarge the definition of X-waves, i.e. Eq. 
©, by introducing the equation: 



(9) 



with K a real constant (k 7^ in the following). Letting 
e = ^{z)'ip{r,t) + i^{z)*'i(j{r,t)* in and setting to zero 
the coefficients of "0 and "0* the following linear system 
is obtained: 



idzfi + K^i + a§x(M + 1^) = 
-idzV + + alxitJ' + 1^) — 0- 



If (/i, ly) — {p, pi) exp(7z) I have 



17 + K + xaf, 
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(11) 



The solvability condition yelds the allowed values for 
the gain 7 = ±-\/— k(k -|- 2xaD. The perturbation 
grows along z if the following inequality is satisfied: 
— k(k -t- Zxoq) > 0. Thus, for a focusing (de-focusing) 
medium, generalized X-waves with —2a^ < k < 
(0 < K < 2aQ) are exponentially amplified. The gain 
Vs K, 7 = \/—k{k + 2x00), plotted in figure 12 has a 
maximum value which corresponds to the most explod- 
ing self-localized packet. 

Now the question arises as to whether or not Eq. Q 
admits solutions resembling X-waves. Closed forms can 
be found proceeding as before: In term of v Eq.Q be- 
comes: 



didittp + Avdvv'4' = Kip, 
whose general solution is 

exp(-\/Kw) exp(\/Ku) 



C2 



(12) 



(13) 
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FIG. 2: Gain 7 Vs k, the parameter identifying the gener- 
alized X-wave, for focusing (thick line) and defocusing (thin 
line) media (ao = 1). 
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FIG. 4: 3D plot of the generalized X-wave ^„ when k 
-A = 1. 




FIG. 3: 3D plot of the generalized X-wave ipK when k = A = 
1. 

A real valued, localised, generalised X-wave is given by 
^^^^( exp(-V^) 

H14|l depends on two parameters A and k; while the first 
determines the decay constant as going far from the ori- 
gin in the (r, t) plane, the latter completely changes the 
shape of the wave. Examples for A = 1 and k, = 1 are 
shown in figure |21 and in figure 0] for A = 1 and k — —1. 
Note that they are similar to the Bessel pulse beams de- 
scribed in . To clarify the differences, I observe that 
the spatio-temporal spectrum develops around the curve 
= fc^ + K, with Lu the angular frequency correspond- 
ing to t and fcj_ the transversal wavenuniber. In figure|Sll 
show the spectra for different k, with k = correspond- 
ing to the simplest X-wave ij^x- The appearance of this 
lines in the spatio-temporal far field is a clear signature 
of the X-waves instability, and can be directly observed 
in the experiments. Note that the spectra resemble the 
shape of the wave in the physical space as a consequence 
of their propagation invariance. 

To show that X-wave instability can actually be ob- 
served in the experiments I numerically solved Eq. 
(0). I considered a gaussian input beam, A — 
Aoexp(-i?V(2Wo2) - TV(2r2)), whose intensity pro- 
file FWHM spot and duration are 70^m and 100/s. 
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FIG. 5: Spatio-temporal spectra (transversal wavevector Vs 
angular spectrum) for different generalized X-waves (thin line, 
K — —2; thick line, k = —1; squares, k = 1; circles k = 4). 
The dashed line is the spectrum of the simplest X-wave (k = 
0). 



Eq. (0) is integrated with reference to fused silica, 
with A = 800nm (no = 1.5,n2 = 2.5 x IQ-^o k" = 
360^^^, SI units ), peak power P = l.5Pc, being Pc — 
(0.61A)27r/(8no"-2) = 2.6MW the critical power for self- 
focusing. In figure El I show the spatio-temporal profile 
and the spectrum (in log scale) after 3 diffraction lengths 
Ldf — ttuoWq/X. Clearly an X-like profile is formed and 
the spectrum shows the features in Fig. [S] 




FIG. 6: Results of the numerical integration of Eq. Q after 
three diffraction lengths ( (Left) Level plots of 101og(|up); 
(Right) Level plots of the square-modulus of the spectrum in 
log scale of u. 
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In conclusion I have shown that a plane wave in Kerr 
media gives rise to linear and exponential amplification of 
X-waves, thus leading to a significat beam reshaping. A 
new class of X-waves is involved in this nonlinear process. 
The reported analysis provides insights for the interpreta- 
tion of the pulse-splitting of focused femtosecond beams, 
and related phenomena, in the same way as MI is relevant 
for solitons generation. Indeed, the spostaneous forma- 
tion of an X-wave can be another explanation of the out 
of axis energy redistribution typically observed. X-wave 
instability can also be an effective approach for the con- 
trolled generation of non-dispersive and non-diffractive 
pulses. 

These results appear to be susceptible of several gener- 
alizations, as considering quadratric nonlinearity or vec- 
torial effects, and have implications in all the fields en- 
compassing nonlinear wave propagation, as acustics, hy- 
drodynamics and plasma physics. For example, X-wave 
instability (as well as nonlinear X-waves) is also expected 
in periodical Bose-Einstein Condensates where Eq. (0) 
holds, being t the direction of periodicity, in the presence 
of negative effective mass. j43, 
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